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ABSTRACT 

A suite of high-resolution N-body simulations of the Magellanic Clouds - Milky Way 
system are presented and compared directly with newly available data from the Hi 
Parkes All-Sky Survey (HIPASS). We show that the interaction between Small and 
Large Magellanic Clouds results in both a spatial and kinematical bifurcation of both 
the Stream and the Leading Arm. The spatial bifurcation of the Stream is readily 
apparent in the HIPASS data, and the kinematical bifurcation is also tentatively iden- 
tified. This bifurcation provides strong support for the tidal disruption origin for the 
Magellanic Stream. A fiducial model for the Magellanic Clouds is presented upon com- 
pletion of an extensive parameter survey of the potential orbital configurations of the 
Magellanic Clouds and the viable initial boundary conditions for the disc of the Small 
Magellanic Cloud. The impact of the choice of these critical parameters upon the final 
configurations of the Stream and Leading Arm is detailed. 
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1 INTRODUCTION 

The progressive collapse and merging associated with hierar- 
chical clustering within a cold dark matter cosmology, while 
do minated by activity at early epochs (redshifts z ~ 2 - 5; 
e.g. iMurali et al1l2002l) . continues to the present-day and is 
readily observable even in the local Universe. In our own 
Milky Way (MW), ongoing satellite dis ruption and accre - 
tion events include the S agittarius dwarf dlbata et aljll994i) , 
the Canis Major dwarf jMartin et alJl2005l and references 
therein), and perhaps the most spectacular of all, the Large 
(LMC) and Small (SMC) Magellanic Clouds (MCs). The 
disruption and accretion of the MCs is perhaps best ap- 
preciated through the nearly circum-Galactic polar ring of 
gas - the M agellanic Stream - emanating fr om the Clouds 
jMathewson et aljll97l iPutman et al.lll99Slh 

Two primary, competing, scenarios have been pos- 
tulated to explain the origin of the Magellanic Stream 
(MS): (i) ram pressure stripping of LMC and SMC gas 
due to the motion of the Clouds through the tenuous 
coronal gas in the Gal actic halo iMoore fc Davis! 11994 
iMastropietro et "all 120051) . This "drag" scenario faces dif- 
ficulty in explaining the Lead ing Arm Feature (LAF) ob- 
served bvJ^taxa^^^L|(|W9^). (ii) tidal disruption of the 
SMC dMurai fc FuiimotolTl980l) . Because of its ability to 
simultaneously produce both trailing and leading streams 
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of gas, the "tidal" m odel has gathered conside r able sup- 
port in the literature ([ Gardin er fc Noguchilll99 li|: ICar dinerl 
| l999tlYoshizawa fc Noguchil200&i hereafter, GN96, G99jand 
lYNO^: respectively). 

Constraints upon theoretical models of the Magellanic 
Stream have improved dramatically wi th the release of the 
Hi Parkes All-Sky Survey (HIPASS: Earnes et alJ 1200 ll) 
dataset. The extent and fine-scale structure of the MS 
can now be appreciated to a level not previously possi- 
ble, in particular the unequivocal e vidence for the exis - 
tence of the Leading Arm Feature iPutman et aljll998lh 
Motivated by HIPASS, we have initiated a program of 
high-resolution N-body modelling of the Magellanic Stream 
aimed at de-constructing the temporal evolution of the 
Magellanic Clouds - Milky Way interaction. Here, we use 
~ 30 times high er reso lution than previous models such as 
lGN96t lG99l and lYNoj 1 . giving us a Hi mass resolution of 
~ 5600 Mq and a Hi flux resolution of 7Jy kms -1 . In ad- 
dition, we construct a detailed Hi map, to compare with 



1 iBekki fc Chibal i2005t) have presented high-resolution simula- 
tions of the Magellanic System, focusing exclusively on the inter- 
nal dynamics and star formation history of the LMC by fixing 
the potentials of both the Milky Way and SMC - i.e. predictions 
concerning the formation and evolution of the Magellanic Stream 
and Leading Arm were not features of their work. The mass res- 
olution (softening length) employed in ou r work is approx i matel y 
a factor of ten (two) greater than that of lBekki fc Chibal |2005|| . 
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the HIPASS data directly and quantitatively. The combina- 
tion of a higher resolution simulation and the simulated H I 
map enables us to argue about more detailed features in the 
MS and the LAF. As a result, this paper shows more con- 
vincingly that the observed LAF and MS can be produced 
by the leading and trailing streams of the SMC, induced 
by tidal interactions with the MW and LMC. Our fiducial 
model is found after an extensive parameter survey of the 
potential orbital configurations of the MCs and the different 
initial condition of the SMC. Based on the parameter sur- 
vey, we also demonstrate how the final features of the MS 
and the LAF are sensitive to the initial configuration of the 
SMC. Our new high-resolution simulations also reveal that 
the tidal interactions create spatial and kinematical bifurca- 
tion in the MS and LAF. We present this prediction based 
on our simulations, however we also demonstrate that the 
existing HIPASS data show such bifurcations. 

The work described here (Paper I) is the first in a series 
of papers aimed at providing a definitive model for the MS 
and the LAF. Here, we concentrate solely upon the effects of 
gravity, describing the SMC disc with collisionless particles, 
and ignore the baryon physics (including hydrodynamics, 
star formation, energy feedback, and chemical enrichment). 
In Paper II, the effects of baryonic physics will be detailed. 

The layout of this paper is as follows: in Section |3|wc 
provide a description of the suite of simulations generated to 
date. In Section |3J we show the results from our best model, 
and compare them with the observational data. In Section[I] 
we demonstrate how the final features of the MS and LAF 
depend upon the orbits of the MCs and the initial properties 
of the SMC disc. Finally, in Section the discussion and 
future directions for our work are presented. 



2 NUMERICAL SIMULATION 

The framework upon which o ur s imulat ions are based paral- 
lels that described bv lGN96l and lYN03t The MW and LMC 
are taken to be fixed potentials, while the SMC is treated 
as an ensemble of self-gravitating particles, in recognition 
of the fact that the MS is thought to originate from the 
tidal disr uption of the SMC disc (e.g. iGardiner et alj| l994l 
hereafte r ICSF94 lGN9rj iMaddison et alJ 120021: lYN03t but 
see also iMastropietro et al.ll2005fl . Theorbits of the MW 
and LMC with respect to the SMC are pre-calcula ted, the 
procedure for which is outlined in Section 12.11 and IGSF94L 
The initial boundary conditions for the "live" SMC model is 
described in Section [2.21 and its evolution explored in Sec- 
tion 



2.1 The orbits of the LMC and SMC 

We assume a spherically symmetric potential for the MW 
halo. The MCs sample a volume of the halo where the gravi- 
tational potential is insensitive to both the expected central 
cusp and the outer regions where the density profile may be 
steeper than an isothermal profile. We hence assume a con- 
stant rotational velocity of V c = 220kms _1 within the MW 
halo. Thus, the potential is described by 



4>a{r) = -V c lnr 

with a mass enclosed within r kpc of 



(1) 



Table 1. Present-day orbital parameters for the SMC and LMC 
from the literature and values we have adopted in this work. 



Parameter 



LMC 
Literature This 
work 



SMC 
Literature This 
work 



"t,GSR E 

f 

b 9 



84±7 2 
281 ±41 EI 
280.46° 8 
-32.89° IH1 
49.43 11 



80.1 
287 

280.46° 
-32.89° 
49.43 



~ 7 3 < 4 
200 ± 100 6 
302.79° El 
-44.30° IH1 
57.oJD] 



7.1 

255 

302.79° 
-44.30° 
57.02 



1 Radial velocity in the Galactic Standard of Rest frame 

j kms" 1 ) 

van der Marel et aTl 12C 



Hard y et al .1 <1989l) 
GSF94I 



° Tangential velocity in the Galactic Standard of Rest 
fr ame fkms" 1 ) 

6 iLin et alJ Il995ll 

7 Galactic lon gitude 

8 iTullvl lll988h 

9 Galactic latitude 

10 Distance (kpc) 

11 iFeast fc Walkeil <1987li 



M 



©■ 



(2) 



M G (< r) = 5.6 x 10 11 f ^ V ( —!—) 

v ; V220kms- 1 / \50 kpc/ 

We assume that there is little disc contribution to the po- 
tential at the typical Galactocentric distances of the LMC 
and SMC (~ 50 kpc). The halos of both the MW and LMC 
are assumed to be invariant for the duration of the simula- 
tion (2.5 Gyr - i.e., a lookback time corresponding to redshift 
z ~ 0.2). Plummer potentials are adopted for both the LMC 
and SMC - i.e., 



6 L , s (r) = GM L , S / [(r - r L , s ) 2 + Kl tS ] 



1/2 



(3) 



where tl,s are the positions of the Clouds relative to the 
MW centre, and Ax,s are the core radii, set to 3 and 2 kpc 
for the LMC and SMC, respectively. In the fiducial model, 
we assume a constant mass of Ml = 2 x 10 10 Mq for the 
LMC, and Af s = 3 x 10 9 Mq for the SMC. 

Adopting these parameters, we next backwards inte- 
grate the orbits of the SMC and LMC, using as bound- 
ary conditions the current tangential and radial velocities 
and positions of the Clouds listed in Tab le ^ (see also 
Murai fc Fuiimotdll980t IGSF94 IGN96L ILin et alJll995l and 
Bckki ct al. 200^). The v alues w e adopted for our work were 
chosen to match those of lGN96l . and are consistent with the 
extant literature. 

The equations of motion of the Clouds about the sta- 
tionary MW are 

d 

rL = Tj— [0s (|r L - r s |) + <f} G (|r L |)] + F L /M L (4) 
and 

f s = -S- [0l (|r s - r L |) + 4>g (|rs|)] + F s /M s , (5) 
Crs 

where the potentials 0l, 4>s and <f>a refer to the Large 
and Small Magellanic Cl ouds, and the Galaxy, respectively 
jMurai fc FuiimotJl 980V Since we do not model the MW as 
live particles, dynamical friction is modelled as per IGSF94L 
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Fl and Fs are the dampening force between the Galaxy and 
each of the Clouds: 

-GM U s 



Fl.s = 



•In (A) 



|r L ,s|~ l rL .sl 

ex P(~!' 2 ) < *!'~ cx P( _:r2 ) :1: ] 



>=[[*L.B| a /V.] : 



(6) 



where In (A) is the Coulomb logarithm. 

Although equation © is a simple analytical formulation 
for dynamical friction that fails to model accurately the or- 
bit of t he Clouds for more than 5 Gyr, according to rec ent 
studies (iHashimoto et alJl2003l: lJust fc Penarrubiall2005l) . it 
gives a good approximation to the predictions of full N-body 
simulations over the ~ 2.5 Gyr period we focus on in this pa- 
per. We note also, the strong effect the C oulomb logarithm 
has on the orbits. lHashimoto et al l j2003f) suggest t h at the 
In (A) value of 3 advocated bv iBrnnevfc^ Tremaind jl987l) 
for the LMC is too large, causing the system to evolve too 
fast. Thus, for this paper, we set ln(A) = 1 to enable direct 
comparison with the previous works. 

Our fiducial model is initiated at a lookback time of 
2.5 Gyrs when the Clouds are at apo-Galacticon (with the 
tidal forces between all three bodies being minimised). The 
subsequent orbit appears in Fig.Q and the distance between 
the Clouds and Galaxy is shown in Fig. [5] 



2.2 SMC model 

The initial dynamical confi guration for the SMC was con- 
structed using GalactICs llKuiiken fc Dubinskilll995l) . We 
employed a bulgeless equilibrium model with a truncated 
exponential disc, and a dark matter King-profile halo. The 
relevant GalactICs parameters for the fiducial SMC model 
appear in Table [5] (a detailed explanation of the individ - 
ual parameters is provided by iKuiiken fc Dubinskil ll995T) . 
The disc possesses an exponential profile with scale radius 
Rd = 3.5 kpc, smoothly truncated beyond R t = 7kpc (with 
95 per cent of both the disc and halo masses being within 
r d,95 ~ ?"h,9S ~ 7 kpc), to give a disc with total radial ex- 
tent (where the face-on surface density of the disc reaches 
1 Mq pc -2 ) of 7.5 kpc (compared to the halo with radial ex- 
tent ~ 14 kpc). The rotation curve peaks at ~ 2 kpc with a 
velocity of ~ 45kms _1 , and turns over to become approx- 
imately constant, giving a total SMC mass of 3 x 10 9 Mq, 
with the disc mass being 1.5 x 10 9 Mq. The central velocity 
dispersion of the disc was chosen to be near the current H I 
velocity dispersion of the observed SMC (Fig. |SJ, and the 
dark matter halo velocity dispersion is similar to the value 
of ~ 25kms _1 obtained from observa tions of the stellar 
halo carbon stars and planetary nebulae jDopita et aljl985t 
lHardv et alJll98St lHatzidimitriou et al.lll997T) . The Toomre 
Q-parameter at the disc half-mass radius is Q = 1.4. 

Our fiducial SMC model assumes a somewhat different 
scale length and total ex tent for the SMC disc when com- 
pared with earlier studies (IGSF94 IGN961 lYN03l) . We should 
stress that one cannot simply take as initial conditions, ob- 
served parameters of the SMC at the present day, since 
the SMC has (obviously) been significantly disturbed by in- 
teractions with the MW and LMC; we instead use initial 
parameters consistent with other relatively isolated dwarf 
disc galaxies (and integrate forward to ensure the final char- 



acteristics is consistent with that observed today). Dwarfs 
with Hi mass 1.3 - 2.0 x 10 9 Mq (comparable to the SMC) 
have Hi disc scale lengths ranging from 1.6 to 4.4 kpc, and 
smaller optical (i.e. stellar) scale lengths of 0.9 to 2.7 kpc 
iS waters et al.l f2002). These same galaxies have H I discs 
whose face-on corrected Hi density drops to 1 Mq pc -2 at 
a radius between 10 to 12 kpc. 

The reasons for the larger truncation radius compared 
to previous work are twofold. First, we suggest the tidal 
radius of the SMC disc is larger than the 5 kpc proposed by 
GN961 Several Gyrs ago, the MW was likely to have been 
slightly less massive (secular halo growth), the SMC slightly 
more massive (less tidal disruption), and the perigalactic 
distance of the SMC was larger. In the estimated orbit of the 
SMC, at a lookback time of 2.5 Gyr (Fig.|2J, the perigalactic 
distance r p was 60 kpc, and the apogalactic distance r a was 
120 kpc. However, the MW mass enclosed within the larger 
S MC orbit has then i ncreased. The tidal radius of the SMC 
is (iFaber fc Lmlll983» 



M s 



(3 + e)Af G (r p ) 



1/3 



(7) 



where the eccentricity e = (1 — fp/r 2 ) 0,5 , and thus r t = 
6.3 kpc when r p = 60 kpc. We also suggest that it is not un- 
reasonable that the disc was initially somewhat larger than 
the tidal radius when the SMC started to experience tidal 
stripping. 

We performed stability tests on all initial SMC models, 
in the absence of any external potential (the equilibrium 
run), to ensure the initial models were indeed in equilib- 
rium. The SMC models were evolv ed using the GCD+ paral- 
lel tree N-body code described bv lKawata fc Gibson! J2003f) 
for 2.5 Gyr (the interaction run). We encountered only a 
minimal degree of disc heating and newly-introduced spiral 
structure (generally with two symmetric arms). 



2.3 Interaction simulations 

We adopt an SMC-centric non-inertial (but non-rotating) 
coordinate system for our interaction simulations, one in 
which the SMC disc lies on the x - y plane. The orbits of the 
LMC and MW from Section f2.1l are translated and rotated 
to this coordinate system, and the SMC model constructed 
in Section 12.21 evolved within this new coordinate system 
under the influence of the now "orbiting" MW and LMC. 
There are two degrees of freedom for the current inclination 
angle of the SMC, both cu rrently unknown. We define 9 and 
<f> following fig. 1 of lGN96t survey the full range of both, and 
and find that the final features of the simulated MS are sen- 
sitive to this angle. We adopt in the fiducial model an angle 
of (9, cf>) = (45°, 210°), as discussed in Sectional This choice 
leads to a trailing tidal stream with an orientation consistent 
with the observed MS, and a leading arm with shape qual- 
itatively similar to the observe d LAF . This angle is mildly 
different from that adopted in lGN96i. (9,<j>) = (45°, 230°), 
but we found that the new value leads to a better match 
to the HIPASS dataset (data which were not available to 
lGN9d) 

The particles in the SMC are assumed to be collisionless 
and their dynamical evolution calculated with GCD +. The 
acceleration applied to the i-th particle is described JGN96I) 
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Figure 1. The three dimensional orbit of the SMC (black line) and the LMC (blue dashed line) around the Galaxy. The rotation axis 
of the MW is assumed to be the z-axis, and the disc plane is centred on z = 0. The orbits, which are derived by backwards integrating 
to T = — 2.5Gyr, are plotted with arrows indicating each 0.5 Gyr between T = and — 2.5Gyr. The green crosses denote the current 
position of the SMC and LMC. The Solar radius at 8.5 kpc is also shown, and the current position of the Sun at T = is marked by the 
open square. 



Table 2. Fiducial SMC model parameters (sec Kuiikcn & Dubinski 1995 for details of parameters). 

Disc Halo 
Af d Rd Rt z A SR ou t ctr.,0 *o cr q C Ra 

(M ) (kpc) (kpc) (kpc) (kpc) (kms" 1 ) (km 2 s~ 2 ) (kpc) 

2.47 x 10 9 3.5 7.0 0.35 1.1 35 - 5.72 x TO 3 29 1.0 0.1 0.96 
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Figure 2. Distance (upper) and approximate tidal forces (lower; 
on a logarithmic scale in arbitrary units) between the three inter- 
acting bodies in our fiducial model. 
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-G 



rrij (n 



,3/2 



i?N {\^-rj\ 2 + e 2 Y 
+ Fmw (r, — tmw) + Flmc fa — tlmc) 
— Fmw (— rMw) — Flmc ( flmc) > 



(8) 



where e is the softening length and rrij is the mass of the j-th 
particle. The position of r is measured in the SMC-centric 
coordinate frame. The first term is the self-gravity of the 
SMC particles; the second and third terms are the forces on 
the particle resulting from the Galaxy and LMC, and can be 
derived from the respective potentials in equations Q and 



V 2 

F M w(r) = -r-Tjr, 
r 



F L Mc(r) = 



GM L r 



(r 2 +Ktf /2 ' 



(9) 



(10) 



The final two terms arise from needing to correct for the 
integration of the equations of motion in a non-inertial ref- 
erence frame centred on the SMC. 

We use 200,000 disc and 200,000 halo particles to de- 
scribe the SMC. This corresponds to a resolu tion ~ 30 
times greater than that employed bv lGN9fjllG99l and lYNOol 
Such high-resolution allows us to examine features of the 
MS, LAF, and SMC, in a manner not previously possible, 
since smaller fractional differences in particle density be- 
come statistically significant. We adopt approximately equal 
masses for the halo and disc particles of 7.6 x 10 3 Mq and 
7.4 x 10 3 Mq respectively, and employ softening lengths ed,h 
of ~ 65 pc. 
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Table 3. Empirical and simulated characteristics of the SMC. 



Parameter 


Fiducial model 


Surveyed range 


Observation 


References 


Hi Mass within SMC (Mq) 








0.35 


- 0.56 


X 10 9 


(1),(2),(3),(4) 


Hi Mass within MS (Mq) 








n 1 

U.l — 


0.5 x 


10 9 


(5), (6) 


Hi Mass within LAF (Mq) 


_ 






n HQ 
U.Uo 


- 0.06 


x 10 9 


(7), (8) 


Hi Mass within ICR (M m ) 


_ 






U.UO 


- 0.65 


x 10 9 


(9), (10), (11) 


Initial Hi mass (Mq) 


1.5 x 10 9 


1.1 


- 2.0 x 10 9 


0.5 - 


1.8 x 


10 9a 




Stellar Mass within SMC (Mq) 


- 






0.58 


- 1.8 x 10 9 


(12),(13) 


Total Mass within SMC (Mq) 


- 






0.9 - 


2.4 X 


10 9 


(14),(15),(16),(17) 


Initial total mass (Mq) 


3 x 10 9 


3.0 


- 3.6 x 10 9 


1.1 - 


3.6 x 


10 9b 




Initial radius of Hi disc (ra 95; kpc) 


7 


4 - 


7 


10 - 


12 c 




(18) 


Initial scale length of H 1 disc (kpc) 


0.5 x r dj95 = 3.5 


0.2 


- 0.5 x r d a 5 


1.5 - 


AM 




(19) 


Initial scale height of H 1 disc (kpc) 


0.05 x r d>95 = 0.35 


0.025 - 0.075 x r dj95 










Initial velocity disp. of Hi disc^ (kms -1 ) 


25 


25 


- 35 


~ 25 






(20) 


Initial radius of DM halo (kpc) 


1 X >"d,95 = 7 


1 - 


1-3 x r dj g 5 











a The sum of the current observed Hi masses within the SMC, MS, LAF and ICR. 

^ The sum of the current total mass within the SMC (obtained dynamically) and Hi masses of MS, LAF and ICR. This is likely underes- 
timated because some dark matter may have been stripped from the SMC, as well. 
c The observed values for SMC-like dwarfs (see text for discussion). 
^ Central velocity dispersion. 

Re ferences: (1 ) iBruns et alj <2005l). ( 2) IStanimirovic et alj Il999fl. (3) IStanimirovic et al.U2004l) who includes He mass, 



f4l lPutman et al.U2003l) . (5) |Putman et alj i2003t>. f 6) lBruns et all <2005t) . C7") |Bruns et alj <2005l) . (8) Upper limit is obtained by 
summing LAF (b < 0°) mass i n lPutmari^talJ ll998l) w ith H VC clouds EP and W D (bu t not WE) complex es in 

IWakker fc van Woerdenl il99lf). f9")lMuller et alj 20031). fl0)|Putman et all <2003l). (ll)lBruns et alj 120051) . (12) IStanimirovic et alj 
J2004l). ri3HYN05 ~7l4) iDopita et alj Jl985l).7l5) lHardv et alj Jl989l) . Q6) iHatzidimitriou et alj Jl9971) . (17) IStanimirovic et alj <2004h . 
( 18) IS waters et al I J2002I) . ( 19) IS waters et al.1 J2002I) . C20) this paper (Fig. |SJ 



3 FIDUCIAL MODEL 

In this section, we show the results of our fiducial model, 
the basic parameters for which are listed in Section [5] and 
Table |3 Since we are primarily interested in the formation 
of the MS and LAF, which we assume are both stripped H I 
gas components from the SMC disc, our analysis focuses on 
the particles which were initially within the original SMC 
disc. In what follows, our simulation products are compared 
closely with the empirical HIPASS dataset. 



3.1 Evolution 

In Fig. [3] we present a series of snapshots in the Galactocen- 
tric coordinate system, with the or bit of the MCs overlayed. 
Consistent with earlier models fe.g. lONQfilYNOsT) . when the 
SMC experienced a close encounter with the MW and the 
SMC ~ 1.5 Gyr ago, the edge of the disc of the SMC began to 
be drawn out, which formed the tidal features that later be- 
came the LAF and MS. By T = —1.0 Gyr, at the subsequent 
apo-Galacticon, the LAF becomes more prominent, whilst 
the MS was still under development. By T = —0.3 Gyr, much 
of the initially stripped material in the leading tidal arm had 
been pulled back into the inter-cloud region (ICR), and the 
material still in the LAF was brought within 3 kpc of the 
solar circle. By this time, the MS and LAF morphology re- 
semble that seen today. The next encounter with the MW 
and SMC at T = —0.2 Gyr caused little obvious consequence 
to either the MS or LAF. It did however, cause the disper- 
sion of much of the material that had been within the SMC. 
Much of this material either ended up in the ICR (although 
most of the ICR material was already in an ICR structure 



before this event), or contributed to the large spread in ra- 
dial extent of the SMC. At the current time, the ICR extends 
radially from ~ 30 to ~ 80 kpc, and the SMC from ~ 45 to 
~ 60 kpc. 

We also ran the simulation "forward" in time to T = 
+0.25 Gyr and found that the ICR undergoes mass loss 
with material being dragged out into two tidal tails, sep- 
arated radially and kinematically from the main MS and 
LAF features. Thi s is consistent with what is expected by 
iBriins et all (120051) . By T = +0.25 Gyr, the material in the 
ICR has been dragged 55 kpc towards the MS in the plane 
of the sky, giving an angular separation of 35° between the 
SMC and tip of the ICR. 

3.2 H 1 column density distribution 

In this section, we compare the present-day Hi column 
density distribution between the empirical HIPASS dataset 
iPutman et aljll998l) and our simulation results. Fig. ^dis- 
plays the Hi column density map on a Zenith Equal Area 
(ZEA) projection. Here, we arbitrarily define the regions cor- 
responding to the MS, LAF, ICR, and SMC as shown in the 
right panel of Fig. [I] We con verted the ob s erved 21cm flux 
to Hi column density, after iBarnes et al. (2001). To con- 
struct Hi column density maps from the simulation results, 
we assume that the disc particles in the SMC are purely 
gaseous, and the H 1 mass fraction is 0.76. The column den- 
sities within the SMC and ICR region will be somewhat 
overestimated, as we neglect currently any associated stellar 
and ionised components. On the other hand, Fig. |S] demon- 
strates that most particles in the MS and the LAF at T = 
originate in the outer edge of the initial SMC disc, where 
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Figure 3. The particle configuration in the Galactocentric coordinate system at different timesteps. Different colour dots represent the 
particles which end up in different components at T = as shown in the bottom-right plot. See the text and Fig. [3] for the definitions of 
the MS, LAF, ICR, and SMC. 



there was likely a lack of stars. To date, searches for a puta- 



tivc stellar component to the MS have 


proven unsuccessful 


(e.s. |PhilirJll976bllal lRecillas-Cruzlll98i 


Brucck & Hawkins 


1983). Regardless, the Hi fraction of 0.76 remains techni- 



cally an upper limit, and thus our predictions should con- 
sidered as upper limits. For these reasons the comparisons 
between the simulation and the HIPASS data focus mainly 
upon the properties of the MS and LAF. 

Fig.0sh.ows that in our fiducial model the gross features 
of the observed MS are reproduced, and the LAF appears as 
a consequence of tidal int eraction s. This confirms previous 
studies, such as GN96 and IYNOoI which suggested that the 
MS and LAF features originate from gas stripped from the 
SMC disc by the tidal interaction with the LMC and MW. 

The left panel of Fig. |7| itself a new representation of 



the HIPASS dataset, demonstrates that the observed leading 
arm extends above the Galactic plane to latitude b ~ 30°. 
While the full extent of the LAF above th e plane remains 
a matter of debate (e.s. iBriins et al]l2005l) , the metallicity 
of its gas is consist ent with an SMC origin (iLu et alJll998l : 
iGibson et all200Cft . supporting the tidal disruption scenario 
for the MS (and LAF). Furthermore, the observed Hi cloud 
distribution seems to show "a kink" near (I, b) — (310°, 0°), 
where two further components of the LAF are seen north of 
the G alactic plane (labelled LAF II and III in lBriins et alJ 
2005). Although the exact position of the kink is inconsis- 
tent with the empirical data, our simulation does naturally 
predict its existence. 

The tail of the observed MS shows spatial bifurcation 
near (l,b) = (300°, -70°) and (l,b) = (80°, -55°) in Fig-H 
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Figure 4. Hi column density map (i.e. zeroth moment map) of the empirical data (left) and our fiducial model (right). A full sky Zenith 
Equal Area projection centred on the South Galactic Pole is applied, to preserve flux. Blue lines delineate regions used to quantify masses 
in Section lij.21 In the right panel, the enclosed regions with labels of MS, LAF, ICR, and SMC are arbitrarily chosen, and are used in 
Fig. |3 The current positions of the SMC and LMC are represented by black and blue ellipses, and in the right panel, the past and future 
1 Gyr histories of the SMC and LMC orbit are denoted by black (solid) and blue (dashed) lines, respectively (the future 1 Gyr is shaded 
in a lighter colour). Lines of constant Magellanic Longitude are drawn in black, with the LMC lying on the 0° meridian. 




-15 -10 -5 5 10 15 -15 -10 -5 5 10 15 

x smc ( k P G ) *smc ( k P c ) 

Figure 5. The distribution of the particles which end up in the MS (left) and LAF (right) at T = in the face-on view of the initial 
(T = -2.5 Gyr) SMC disc. For the definition of the MS and LAF see Fig. El 



with the two component s forming an apparen t twisting dou- 
ble helix-like structure jPutman et al.l Eo03K T his bifu rca- 
tion is not apparent in previous studies such as GN96 and 
lYN03t higher resolution simulations enable us to study such 
subtle features. Since this might be further evidence of the 
tidal interaction between the LMC and SMC, we return to 
this issue later. 

An advantage of the present work is that the H I data 



available to co nstrain th e mode ls are significantly improved 
beyond that of lGNQtl or lYNOfl As in those previous stud- 
ies, while the gross features of the observed MS and LAF 
are reproduced by our fiducial model, there remain subtle 
discrepancies between simulations and data. While the sim- 
ulated MS is both broader and more extended than that 
observed (and hence the mean column density is somewhat 
lower than that encountered in the HIPASS dataset), the 
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Figure 6. The heliocentric distance, as a function of the Magel- 
lanic longitude, derived from our fiducial model. The black (solid) 
and blue (dashed) lines show the past and future 1 Gyr histories 
(the future 1 Gyr is shaded a lighter colour) of distance and Mag- 
ellanic longitude for the SMC and LMC, respectively. 

derived Hi gas mass (assuming a heliocentric distance of 
57kpc) from HIPASS is within 50 per cent of that of our 
simulated fiducial model (3.5 x 10 s Mq and 2.4 x 10 8 M , 
respectively). Fig. [(J shows the heliocentric distance of the 
simulated MS and LAF against the Magellanic longitude, 
#mc (as defined in Wakkcr 200ll: lines of constant #mc are 
shown in Fig. . This demonstrates that the distance to the 
simulated MS is not constant, instead increasing across its 
length. Hence negating the distance ambiguity by obtain- 
ing a total flux gives us a fairer comparison than a total 
mass. Doing so results in an observed MS total Hi flux of 
4.8 x 10 5 Jy kms -1 , while the fiducial simulated MS has a 
total Hi flux of 2.3 x 10 5 Jy kms -1 , still within a factor of 
~ 2 of the HIPASS data. 

Conversely, the simulated LAF has a predicted H I mass 
and associated flux (7.3 x 10 7 Mq and 7.9 x 10 4 Jy kms -1 , 
respectively) both factors of ~2 greater than that in- 
ferred from the HIPASS dataset (3.5 x 10 7 Mq and 4.6 x 
10 4 Jy kms -1 , respectively) 2 . One clear difference between 
observation and simulation is that of the geometry of the 
LAF, in particular that of the projected deflection angle be- 
tween the LAF and a Great Circle aligned with the MS 
proper (Fig. 0. In addition, the simulated LAF extends 
above the Galactic Plane beyond that observed. We will 
discuss possible solutions to these apparent problems in Sec- 
tion 

3.3 Hi kinematics 

Fig. |H1 shows the first and second moment maps for both 
our fiducial simulation and that derived from the HIPASS 

2 We have defined the LAF and MS regions differently between 
the empirical and simulated datasets, to account for the geomet- 
rical differences, such as angle, width, and length, between the 
two. 



dataset. In order to remove the large velocity gradient along 
the MS (which acts to obscure fine kinematical details within 
the maps), the first moment map is shown as the distribution 
of the velocity of «gsr — «ms, where vms is the mean trend of 
velocity across the observed MS and LAF, and is defined in 
terms of two fitted Fourier components, vms = 86 sin(#MC + 
2.5°) — 92sin(2(#MC — 26.3°)). The second moment map is 
the distribution of the velocity dispersion of H I gas. 

In the first moment map, the observed velocity trend 
along our simulated MS is shown to be globally consistent 
with that of the empirical data. Fig. [5] displays i>lsr against 
the Magellanic longitude, 9mc, and demonstrates perhaps 
more clearly that the mean velocity of the simulated MS is 
consistent with that observed. In contrast, the line-of-sight 
velocity of the simulated LAF is significantly larger than 
that observed, although it does follow the general trend of 
decreasing velocity at 6mc > 0. 

The second moment maps seen in the lower panels in 
Fig. |H| indicate that the velocity dispersion of the simulated 
LAF is roughly consistent with, although slightly higher, 
than that observed. However, the velocity dispersion of the 
simulated MS is somewhat greater than that inferred from 
the HIPASS dataset. This might be due to the neglect of 
gas dissipation via radiative cooling in the current suite of 
simulations, as will be examined further in Paper II. 

We end by drawing attention to the evidence of a bi- 
furcation in wlsr (right panel of Fig. |UJ within the MS and 
LAF. The next section discusses this bifurcation in more 
detail. 

3.4 Spatial and velocity bifurcation 

Our fiducial model shows a bifurcation in the MS, as seen 
in the ZEA spatial distribution of the H I column density of 
Fig. 2] Within the simulation, this bifurcation occurs both 
radially and tangentially. Fig. 01 displays the tangential bi- 
furcation, and shows that there are two stream components 
that appear to follow a twisting topology governed by the or- 
bits of the MCs - both the orbits and the bifurcation seem 
to cross at (l,b) ~ (45°,— 80°). Radially, the head of the 
Stream, just behind the MCs, has two components (Figs. [3] 
and|nj. At the tail of the MS, there is a third component, well 
separated from the other two streams — extending from 170 
to 220 kpc. It is visible only through its heliocentric distance 
from the Sun, and coincides with the (I, b) position of the tip 
of the main streams. Also apparent in Figs. HJ and |7| are the 
radial and tangential bifurcations of the LAF, respectively. 

The right panel of Fig. demonstrates that the bifur- 
cation of the MS appears also in «lsr- Plotted in Fig. |2l 
are lines showing the history (and future) of the MCs 
wlsr- Only one of the bifurcated components follows the 
orbit of the Clouds (primarily the LMC), while the other 
possesses a higher velocity. Interestingly, there is a second 
velocity component at the position of the SMC (as well 
as at the h ead of the MS and the L AF). We remind the 
reader that iMathewson fc For 1 (11984) observed two veloc- 
ity components within the SMC itself, with a separation of 
~ 50 kms -1 . While the separation of our two velocity com- 
ponents is much larger, it might indicate that the two ob- 
served velocity components are caused by a similar process. 

Snapshots of the simulation, similar to those seen in 
Fig. [3] hint at the origin of these bifurcations. Prior to the 



© 2005 RAS, MNRAS OOO.ITIITIl 



N-body Simulations of the Magellanic Stream 9 




Figure 7. Hi column density map of the LAF region derived from the HIPASS dataset (left) and that of our fiducial model (right) 
with ZEA projection centred on the LAF. The future orbits of the MCs are denoted by a lighter shading, demonstrating the similarity 
between the orbit of both MCs and the orbit of the LAF. 



first major peri-Galacticon at T = — 1.5 Gyr, an encounter 
with the LMC at T = —2.2 Gyr drew the particles from 
the tip of the SMC disc closest to the LMC (and furthest 
from the MW), which resulted in the particles that even- 
tually came to reside in the most distant MS component, 
at a Galactocentric distance of 170 to 220 kpc. The MS 
particles did not become "distinct" from the SMC proper 
until T = -1.5 Gyr; at T = -1.05 Gyr, the MS then re- 
ceived an impulse from an LMC encounter, which caused 
the spatial bifurcation of the MS. The MS was then given 
a "kick" by the LMC at the subsequent apo-Galacticon at 
T = -0.55 Gyr. This encounter at T = -0.55 Gyr resulted 
in the MS being broken into two kinematic components, re- 
sulting in the apparent velocity bifurcation. 

On the other hand, a portion of the LAF comes from 
the particles on the same side of the disc, but closer to the 
SMC centre at the T = —2.2 Gyr encounter. The opposite 
side of the edge of the disc mostly consisted of particles that 
end up at the current time in the ICR. At T = —0.9 Gyr, 
the LMC passed through the LAF, splitting it into two bi- 
furcated radial components. 

Our models seem to create naturally both spatial and 
velocity bifurcations, via tidal interaction between the LMC 
and SMC. If this is the case, any observed bifurcation fea- 
tures would be strong evidence supporting the tidal forma- 
tion scenario for the MS. Observationally, the left panel of 
Fig.|l|shows that there is spatial bifurcation in the observed 
Hi distribution. However, from this data alone, it is difficult 
to ascertain from the Wlsr distribution of Fig. [5] whether 
there is an observed veloc ity bifurcation. Nev ertheless, both 
Fig. |3] of this work, and iBriins et alJ i2005l) . with data of 
higher velocity resolution, show evidence for a bifurcation 
in ulsr along the MS, with the two components observed in 
this work being separated by approximately 100 km s -1 at 
8mc ~ 100° , and two components being visible in fig. 3 of 
IBriins et alJ <l2005l) between the interface region and Galac- 
tic Plane. There is also a hint of bifurcation (with separation 
of ~ 100 kms -1 ) in the LAF at 6 M c ~ 45°. Unfortunately, 
the current quality of the observational data is not enough 
to lead to a firm conclusion. Higher quality data cubes with 
improved velocity resolution will provide critical information 
concerning this apparent bifurcated structure. 



4 PARAMETER DEPENDENCES 

As explained in Section |2] our simulations involve several 
parameters which are not well-constrained by current ob- 
servations. In this section, we briefly demonstrate how the 
final configurations of the MS and LAF are sensitive to 
these free parameters. We varied parameters over a wide 
range of parameter space in our survey, to obtain our fidu- 
cial model. The parameter survey was performed at a lower 
resolution (25,000 disc and 25,000 halo particles) than the 
final fiducial model. We have confirmed that in the fiducial 
model the results of the lower-resolution simulation are con- 
sistent with the higher-resolution simulation, although, as 
presented above, the choice of a higher resolution enables us 
to discuss more detailed features. For example, in Fig. 1101 
we show the H I column density map for the lower-resolution 
model with identical parameters as our fiducial model. The 
distribution of the H I column density is roughly consistent 
with that shown in Fig. 2] although the higher-resolution 
models affords an improved examination of the finer-scale 
structures intrinsic to the simulated streams. The bifurca- 
tion alluded to in the previous section is an example of such 
a feature, and is perhaps best appreciated through an in- 
spection of Fig. 1111 where the flux "excess" of the high- 
resolution model with respect to the lower-resolution model 
is presented in the vlsr vs. 9mc plane, and the bifurcation 
becomes becomes readily apparent. 

Most of the parameter space surveyed is summarised in 
Tabled The scale height of the SMC disc, the ratio of SMC 
disc mass (Hi mass in Table |3J to SMC halo (DM) mass, 
the velocity dispersion of the SMC disc, and the initial and 
final total mass of the SMC (varying the mass of the SMC 
for the purposes of the orbit calculation) were found not 
to be important to the evolution of the system. Both the 
velocity dispersion, and the ratio of H I to DM mass of the 
disc, simply scale the velocity dispersion and quantity of H I 
found in the final Streams and Clouds in a linear fashion. 

In the parameter survey, the orbits of the MW, LMC 
and SMC were derived assuming the mass of the SMC is 
constant. However, we found that the SMC bound mass de- 
creases from 3 x 10 9 to 1.5 x 10 9 Mq approximately linearly 
with time between T = —2.25 Gyr (near the first interaction 
between the LMC and SMC), and T = 0. Such mass- loss 
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Figure 8. First (upper) and second (lower) moment maps derived from the HIPASS dataset (left) and that of our fiducial model (right). 



may affec t the orbit and the final featu res of the MS and 
LAF (e.g. lZhaolEoollKnebe et al.ll2005|) . Thus, we explore 
the effects of decreasing the mass of the SMC linearly with 
time from 3.0 x 10 9 to 1.5 x 10 9 Mq. We found that there is 
little change to the orbits, since the LMC, rather than the 
less massive SMC, primarily determines the orbit of both 
bodies in equations @ and @. As a result, we confirmed 
that the final features of the MS and the LAF are also not 
influenced heavily by the time evolution of the SMC mass. 
Thus, in all other simulations in the parameter survey, we 
use the orbits predicted with no evolution of the SMC mass. 
We now highlight the influence of the most important input 
parameters. 



4.1 The initial scale length of the SMC disc 

We found that the final H I distributions are sensitive to the 
scale length of the initial SMC disc. Fig. ll2l shows the column 
density map of a model with a smaller scale length (1.4 kpc) 
for the initial SMC disc. This reduced scale length results 
in a lower total Hi flux for the MS, as the initially more 
concentrated SMC mass distribution results in less material 
being stripped from the disc. Since the fiducial model does 
not have a total H I flux high enough to match perfectly the 
observations (Section 13.21 , we conclude that reducing the 
SMC scale length is not appropriate. It is also worth noting 
that the smaller scale length model leads to a less significant 
bifurcation in the simulated MS. 
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Figure 9. The line-of-sight velocity with respect to the Local Standard of Rest «LSRi as a function of the Magellanic longitude 6mCi 
derived from the observed HIPASS dataset and that of our fiducial model. The black (solid) and blue (dashed) lines show the past and 
future 1 Gyr histories (the future 1 Gyr is shaded a lighter colour) of i>lsr an d #MC (where the longitude is denned with respect to 
the current position of the LMC) for the SMC and LMC, respectively. A reference line (dotted green) is drawn to show the velocity 
subtracted equally to yield the top panels of Fig. [H] — see text for details. 




the finer features visible. Contours are plotted at 5 X 10 39 and 
Figure 10. Hi column density map for the lower-resolution sim- 5 X 10 40 (cmkms -1 ) -1 . 

ulation with the same parameter set as the fiducial model. 



To quantify the difference between the models, we have 
measured the H I masses and total fluxes within the regions 
delineated as the MS in Fig. [I] As a result, the low- resolution 
simulation of the fiducial model has an H I mass in the MS of 
2.4 x 10 8 M Q and total Hi flux of 2.5 x 10 5 Jy kms -1 . The 
low-resolution simulation has a very similar mass to that 
of the high-resolution case (see Section 13.211 . On the other 
hand, the small scale length model leads to a MS Hi mass 
and total flux of 1.8 x10 s M and 1.9 xlO 5 Jy kms" 1 respec- 
tively, which is significantly smaller than that of the fiducial 



model. Since the fiducial model has a mass which is some- 
what lower than that inferred from the HIPASS dataset, we 
again conclude that the models derived with the reduced ini- 
tial scale length for the SMC perform worse than the fiducial 
selection. 

The H I surve y of late-type dwarf galaxies by 
ISwaters et al.l <l2002l) suggested that the range of scale 
lengths of the gas disc for galaxies which have a similar H I 
mass to that of the SMC is 1.5 - 4.5 kpc. The scale length 
of our fiducial model (3.5 kpc) appears reasonable, and thus 
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Figure 12. Hi column density map for models with a smaller 
scale length (1.4 kpc) for the initial SMC disc. 



the progenitor of the SMC is considered to have had a large 
Hi disc before the SMC fell towards the MW. 

We note that choosing a smaller scale radius is simi- 
lar to setting a small truncation radius, in that the SMC 
material is concentrated more strongly towards the centre 
of the SMC, where it is more difficult to tidally strip. We 
obtain similar results to the above when we choose smaller 
SMC truncation radii, such that the stream is substantially 
retarded when the truncation radius is reduced from 7 kpc 
to 4 kpc. In this situation, the MS Hi mass and total flux 
are reduced to 7.2 x 10 7 M Q and 7.0 x 10 4 Jy kms~\ The 
LAF mass is scaled in a similar manner in both cases. 



4.2 The inclination angle of the SMC 

As mentioned in Section 12.31 the current inclination angle 
of the SMC is unknown, and the initial disc angle 9 and (j> 
are free parameters in our simulations. Fig. 1131 shows the 
Hi column density distribution for the models with (8, <f>) = 
(30°, 210°) and (45°, 230°), to demonstrate how these angles 
affect the final distribution of the MS and LAF. We remind 
the reader that we use (9,(f>) — (45°, 210°) in the fiducial 
model, and therefore small differences of only 20° in the 
initial inclination angle have quite a marked effect on the 
details of the final distribution. If accurate observations of 
the current inclination angle of the SMC were to be made, 
it would provide a strong constraint on any putative model 
of the formation of the MS. 



4.3 The mass of the LMC 

Another unknown para meter is that of the mass of the LMC 
itself. At the time of the lGN9fJ study, the mass of the LMC 
was believed to be ~ 2 x 10 10 Mq (e.g. ISchommer et all 
1992). However, r ecently, some authors claim the lower mass 
of the LMC (e.g. Ivan der Marel et alJl2002l who suggested 



that the mass of the LMC within R < 8.9 kpc is about 9 x 
10 9 Mq). Motivated by such claims, we also ran models with 
varying orbital parameters, and in particular, a lower LMC 
mass. Fig. 1141 shows the result of a model with an LMC 
mass of 1.5 x 10 10 Mq. Even such small differences in the 
LMC mass cause a large change in the orbits of the LMC 
and SMC. Since the MS follows the past orbit of the MCs, 
the angle of the MS is radically different between the model 
with small LMC mass, and both the fiducial model and the 
observed MS. This result suggests that such a low mass LMC 
is unlikely, if the MS is the result of tidal interactions. 



5 SUMMARY AND DISCUSSIONS 

We have carried out high-resolution N-body simulations of 
the history of the SMC disturbed by tidal interactions with 
the MW and LMC. We have surveyed most of the possi- 
ble parameter space for the SMC orbit and the properties 
of the initial SMC disc, and found the best model. The in- 
creased numerical resolution of ~ 7Jy kms -1 per particle 
for the Hi flux, ~ 30 ti mes higher than the previous stud- 
ios made it possible for us to exam- 
ine the detailed features of both leading and trailing tidal 
streams. Taking advantage of this higher resolution, we for 
the first time made a direct and quantitative comparison of 
the simulation results with newly available high quality ob- 
servational data from HIPASS. We convolved the HIPASS 
dataset with the identical software tools used to analyse 
the simulated datasets, and compared the results in iden- 
tical manners. Such comparisons confirm the conclusions of 
previous studies - that the existence of the LAF and MS 
can be explained by the leading and trailing streams of the 
SMC created by tidal interactions with the MW and LMC. 
However, our quantitative comparison have revealed extant 
problems with the models (some minor; some more signifi- 
cant). We found that even in our best model (1) the shape 
of the MS is too extended in both width and length (Sec- 
tion l3.2|l : (2) the total H I flux of the MS is too low, and thus 
the mass of the MS is too low (Section l3.2H ; (3) the angle the 
LAF emanates from the MCs is not entirely consistent with 
observations (Section I3.2fl : (4) and nor is its total Hi flux 
consistent with observations (Section 13. 2H ; (5) the velocity 
dispersion of the MS is too high (Section 13.31 : and (6) the 
line of sight velocity of the LAF is too high (Section 13.3^ . 

These problems suggest that additional physics may be 
required to explain the observed properties of the MS and 
LAF, quantitatively. Obvious physics which this study ex- 
cludes are gas physics, such as hydrodyn amics and dissi- 
pation by radiative cooling. Simulations in lMaddison et ail 
(2002) indicate that dissipation causes the MS to become 
narrower, which might lead to a reduce d velocity dispersion. 
Ram pressure llMoore fc Davisl Il994t iMastropietro et alJ 
2005) is another physical process which our numerical model 
does not currently ta ke into account. Ram pressure (or drag 
- see lGardinerlll999l) is expected to shorten the leading arm 
and increase the gas density in the MS, which should help 
to solve the deficiencies of the model at recreating the LAF 
and MS with their correct shapes and densities. If it were 
to also bring the MS significantly closer than the 57 kpc 
fixed for the observational data within this paper, then the 
mass of the observed stream and total flux of the modelled 
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Figure 13. Hi column density map for the model (0,0) = (30°, 210°) (left) and (45°, 230°) (right). The current positions of the SMC 
and LMC are represented by a black ellipse (the ellipse is the projection of the SMC disc at the given angle) and blue circle respectively. 




Figure 14. H I column density map for the model with an lower 
LMC mass of 1.5 X 10 10 Mq. The past and future 1 Gyr histories 
of the orbit resulting from this differing LMC mass are presented 
by black solid and blue dashed lines for the SMC and LMC re- 
spectively 



modelled the MW potential as spherically symmetric with 
a constant rotational velocity of 220 km s -1 , any deviations 
from this (such as the unknown contribution from the MW 
disc; e.g. iFich fc Tremaineill99lh would affect our modelled 
LAF in particular. We are in the midst of introducing these 
physical processes into our simulations, and will report upon 
their respective effects in Paper II. 

Another benefit of the higher-resolution simulations 
performed here is the identification of the bifurcation of the 
MS and LAF both spatially and kinematically. Our simu- 
lations predict that if the MS is created by a tidal inter- 
action with the LMC, the bifurcation would appear both 
in the H I column density map and the line-of-sight veloc- 
ity versus Magellanic longitude plane. Current observations 
are consistent with the existence of a spatial bifurcation in 
the Hi column density map (left panel of Fig. [IJ. In the 
velocity versus Magellanic longitude plane, it is difficult to 
make a firm conclusion from the current HIPASS data, al- 
though it is interesting that the SMC it self is found to con- 
sist o f at least two velocity components iMathewson fc Fordl 
Il984h . perhaps caused by the same tidal disruption processes 
that form the bifurcation in our models. Observations with 
high velocity resolution and sensitivity may be able to test 
our prediction of bifurcation within the MS. If confirmed, it 
would provide strong evidence that the MS was created by 
tidal interactions. 



stream would be increased, bringing them back into agree- 
ment. However, drag might lead to more problems for the 
length of the MS. Another possibly important mechanism 
is supernovae feedback which might aid in ejecting gas from 
the SMC and/or LMC, and help to increase the gas density 
in the MS and LAF. Finally, we note that the LAF passed 
very close to the centre of the MW 0.2 Gyr ago, and since we 
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